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We review recent results on the potential energy landscape 
(PES) of model liquids. The role of saddle-points in the PES 
in connecting dynamics to statics is investigated, confirm- 
ing that a change between minima-dominated and saddle- 
dominated regions of the PES explored in equilibrium hap- 
pens around the Mode Coupling Temperature. The structure 
of the low-energy saddles in the basins is found to be simple 
and hierarchically organized; the presence of saddles nearby 
in energy to the local minima indicates that, at non-cryogenic 
temperatures, entropic bottlenecks limit the dynamics. 



The study of the properties of the free energy land- 
scape and/or of the potential energy surface (PES) and 
their connection with the slow dynamics in undercooled 
liquids and disordered systems is a topic of current re- 
search [Q. The original idea of Goldstein of consid- 
ering the atomic dynamics of undercooled liquids as the 
incoherent superposition of a vibrational motion confined 
in the basin of attraction of local minima of the PES, in- 
terspersed by jumps among the basin of different minima, 
was later developed in a formal theory by Stillinger and 
Weber These authors define the basin of attraction 
of a given local minimum of the PES as the set of all 
that points in configuration space who are linked to the 
minimum through a steepest descent path. It is obvious 
that - according to this definition - the basins of attrac- 
tion generate a complete partition of the configurational 
space. Consequently, Stillinger and Weber write the par- 
tition function as 

Z = Y,Zb (1) 

{B} 

where {B} indicates the set of all the basins and Zb is the 
partition function restricted to the basin B. The points 
belonging to the borders {dB} of the basins are safely 
excluded from the calculation of the partition function, 
as they form a set of zero measure. 

In this way Stillinger and Weber succeeded to write 
the free energy of a system with a continuous PES as the 
sum of the free energy of a typical basin Fbasin (T, p) and 
a term accounting for the multiplicity of basins accessible 
by the system fl{T, p) — exp(S'con/(2^, p)/kB)- Excluding 
the local minima corresponding to crystalline structures, 
they arrive to an expression for the total free energy of 
the liquid phase: 

Fliquid — Fbasin{T, p) — TSconf{T, p) (2) 



that can be used to predict, for example, the melting tem- 
perature Tm as that temperature where Fnquid becomes 
equal to the free energy of the crystal Fcrystai ■ 

According to the Stillinger and Weber idea, the liq- 
uid thermodynamics can be derived from the knowledge 
of the basins of the PES. On the contrary, the proposal 
of Keyes to use the instantaneous normal mode (INM) 
theory to predict not only the very short time dynam- 
ics but also the long time diffusion processes jij gives 
relevance to the process of crossing the borders of the 
basins and, therefore, to the structure of these borders. 
In the INM approach, the local curvature of the PES is 
calculated along equilibrated molecular dynamics (MD) 
trajectories. The idea is to reconstruct the global shape 
of the landscape explored from local information. For 
example, if the system is near enough to a local mini- 
mum, the local curvatures of the PES will be all posi- 
tive and similar in value to the curvatures of the nearby 
minimum. On the other hand, points near the borders 
between basins will have negative curvatures along some 
directions. In fact, the potential energy has a maximum 
along any direction crossing a border (Fig. 0a). Although 
the exact functional form is system dependent Q, the 
INM approach suggests that the diffusion coefhcient D is 
a monotonic function of the fraction of unstable modes 
<fu>=<Nu>/N, where <iV„> is the number of nega- 
tive eigenvalues of the Hessian matrix evaluated along 
the actual MD trajectory. 

A more quantitative relation between the number of 
negative curvature of the PES and the diffusion coeffi- 
cient stems from the observation that not all the direc- 
tions with negative curvatures lead to a change of basins: 
even in the crystalline state, a non-zero <fu> is found 
due to the anharmonicities of the potential. Depend- 
ing on the different authors and on the different model 
studied, different methods to discriminate the fraction 
<fdiff> of truly diffusive directions from the whole set 
of negative curvatures have been proposed. A non zero 
value of <fdiff> indicates that the trajectory followed 
by the system is confined to regions of the configuration 
space near the borders of the basins |^ . 

The study of the diffusion coefhcient in undercooled 
liquids naturally leads to a comparison with the Mode- 
Coupling Theory (MCT) §. Indeed, the ideal MCT pre- 
dicts the existence of a finite dynamical critical temper- 
atures Tmct at which D = 0. The ideal mode-coupling 
transition can be derived under the restricted assump- 
tions of the gaussianity of the fluctuations [D . The INM 
approach reveals that at Tmct there is a change in the 
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temperature behaviour of the quantity <fdiff>' it be- 
comes numerically zero (Fig. ^3). Therefore one can con- 
clude that at Tmct the system ceases to be near the bor- 
ders of the basins and starts to spend most of the time 
close to the minima. 

On a general ground, one has to understand why the 
system is not already close to the minima at undercooled 
temperatures Tmct <T < Tm, or, in other words, why 
the slowing down of the diffusion processes takes place 
even if the system spend most of the time close to the 
basin borders. The analysis of Kurchan and Laloux ||] on 
the properties of borders in high-dimensional configura- 
tion spaces indicates a simple reason for this behaviour. 
In a A^-dimensional space, (i.e., the configuration space 
of a system with N degrees of freedom) a basin B has a 
volume at least of the order l'^ , where I is a finite length 
and is the typical distance among configurations in the 
same basin. We can imagine then to give a fictions vol- 
ume to the border dB considering "belonging to dB" all 
the points in a shell of width e and volume ^ eNl'^^^ 
near dB. The fraction of border points is therefore eN/l 
indicating that in the thermodynamic limit N 00 
most of the volume of the basins is located in a shell 
of width e 1/N near the borders. The numerically ob- 
served decrease of <fdiff> with the temperature gives 
us further indications on the structure of the borders. At 
high temperatures, the system populates regions of the 
borders that are in common to many basins, while near 
Tmct the system is near regions of the borders in com- 
mon to few basins; therefore the regions of the borders 
in common to few basins have a smaller volume com- 
pared to regions in commons to many basins. This seems 
to indicate a hierarchical organization of the configura- 
tion space p. Indeed, let's assume that we are able to 
study the potential energy U restricted to the borders 
dB of the basins. This restricted PES is an {N — 1)- 
dimensional manifold where all the stationary points |lC|] 
have their order decreased by one: so, the saddles of order 
one become the local minima of the restricted (A^ — 1)- 
dimensional PES and partition this manifold in basins. 
Again, if the system is at high temperatures, the trajec- 
tories will be localized near the (A — 2)-dimensional bor- 
ders of the (A^— l)-dimensional basins of attractions; so, 
one can iterate the procedure until one consider the rele- 
vant (A^ — i)-dimensional manifold in the nearby of which 
equilibrium trajectories are "almost always" within some 
small distance (with the order i less than or equal to 
N<Uff>). 

Following the previous rationale, in a rough way, one 
could write the partition function of the system as 



(3) 



where Mi is the set of points in common to i basins and 
ZMi is the partition function restricted to the (A^ — i) 
dimensional manifold Mi. The local minima of the po- 
tential energy U restricted to Mi are the saddles of or- 
der i in the unrestricted configuration space. Although 



from the strict point of measure the manifolds Mi have 
zero measure, we know from the INM studies that at 
temperatures T > Tmct the system is in the nearby of 
some manifold Mi(^x) with i{T) > 0. In order to give 
an operational definition of such partition of the space 
in basins of attractions of the saddles, one has to give 
a finite volume to the borders. To this aim, two recent 
studies ||Tl| , h ave used the pseudo-potential introduced 
by StiUingeFpl]: 



$ = |VC/| 



(4) 



It is worth to note that the stationary point of any 
order of the PES are absolute minima of such potential 
(in particular they are characterized by $ = 0), but there 
exist minima of <I> - those with <I> 7^ - that are not 
stationary point of U. This implies that the use of is 
not rigorously correct if the goal is to obtain a partition of 
the configuration space in basins of attractions of saddles 
of different orders. 

The choice of <& provides anyhow a partitioning of the 
configuration space in basins of attraction of its local 
minima. As shown recently i n ||l4| and confirmed by a 
more detailed analysis by us |15[, the majority of local 
minima of are not true saddles but present a 2% frac- 
tion of directions which are indeed inflection points. De- 
spite this, the T-dependence of the number and energy 
location of these "quasi-saddles" carries informations on 
the dynamical slowing down upon cooling There- 
fore, although not all the minima of $ found can be clas- 
sified as "true" saddles we will still assume that 
they are a good sampling of the PES and therefore that 
the "true" saddles maintain the same characteristics (i.e. 
curvatures, number of negative directions, energy) of the 
found local mimima of $; in this respect, we will continue 
using the word "saddles" in this manuscript, keeping in 
mind that we refer to such minima. 

For the case of simple model liquids the order <Ns> 
and the energy of the saddles explored during the MD 
simulations have determined As a result, the 

saddle's order has been found to be a decreasing func- 
tion of the temperature and to become zero at the Tmct 
temperature of the system under considerations. This 
result confirms the role of Tmct as a crossover temper- 
ature from a dynamics near the basins's borders {<Ns> 
> 0) to a dynamics mainly localized inside the basins 
{<Ns> — 0) jl^]. Given this findings, the strong non- 
Arrhenius increase of the relaxation times of the dynam- 
ics approaching T^/ct from above can be explained with 
the decrease of the diffusive direction, whose number is 
proportional to <Ns>. Indeed, although the system does 
not need to cross any energy barrier to "jump the ridge" 
among basin, the exploration of the configuration space 
is slower as the connectivity (i.e. the number of adjacent 
basins that the system can freely reach) among basins 
decreases. 

A direct evidence of situations in which, although the 
energy of the saddles is significatively less then the en- 
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ergy of the system, the system is trapped in the basin of 
attraction of a local minimum, is found in [^ , where the 
connectivity between the crystalline ground state and the 
adjacent local minima is analized for a model of the crys- 
talline compound MgF2- Moreover, it is found in | |T7| ] 
that, at higher energies, saddles connecting several min- 
ima make up a sizable fraction of the transition region. 

A further confirmation of the " entropic" - rather than 
"energetic" - origin of the slowing down of the dynamics 
is given by numerical aging experiments. For Lennard- 
Jones liquids, it has been shown that in the aging follow- 
ing a deep quench, the system follows quasi-equilibrium 
states [ p^ during which short time dynamics is intra- 
basin vibration, while long time dynamics is dominated 
by the slow search of deeper and deeper local minima. 
The work of Angelani et al. shows that, after a short 
transient, the system changes basins crossing saddles of 
order one without any energy activated process. 

The result that at Tmct the value of <Ns> ap- 
proaches zero has lead to the interpretation of the mode 
coupling critical temperature as the threshold below 
which the diffusive dynamics becomes activated with po- 
tential energy barriers that the system needs to over- 
come in order to cross the border between basins [ pO| . 
We show here - after a careful definition of the meaning 
of the term "activated processes" in the framework of 
the PES description of the dynamics - that it is at least 
unnecessary to invoke potential energy barrier crossing, 
and that there is no qualitative change in the way the 
system change basins below or above Tmct- What is 
qualitatively different on crossing Tmct is the vanishing 
probability of occupying a basin of a saddle below Tmct ■ 

Let's first discuss the concept of activated processes 
in the framework of landscape theory. Landscape theory 
considers the motion of a representative point in an TV- 
dimensional configuration space; in the thermodynamic 
limit TV — > oo, the value of thermodynamic quantities 
like the average potential energy U becomes determined 
apart of fluctuations of order l/ViV Although in 

ergodic theory the complexity of the landscape is not 
explicitly considered but as the hypotesized source of 
random couplings enforcing the large-numbers law and 
the ergodicity |22|, the apparatus of the ergodic theory 
approach to the foundation of statistichal mechanics is 
based on the analysis of the motion of the system on iso- 
potential manifolds of the PES. The exploration of such 
manifolds is purely entropic, and no fluctuation in en- 
ergy is considered. Obviously, for N large enough, one 
can always consider a subset of the system with a fi- 
nite number M degrees of freedom. This subset can be 
then considered as a canonical system of M degrees of 
freedom with finite fluctuations of order 1 / VM . In 
general, it is possible to select a subset of coordinates 
{lo, the "reaction" coordinates) relevant to the process 
under consideration, as for example the coordinate de- 
scribing the motion across the basin's border. In this 
case the potential energy of the subset along the reaction 
coordinate goes through a maximum (a saddle point in 



the configurational space) . The other degrees of freedom 
can be considered as a "thermal bath" for the reaction 
coordinate. 

Our aim is now to show that both approaches - the one 
considering the systems as a whole moving on a constant 
potential surface and the one that focus the attention 
on the reaction coordinate and consider the rest of the 
system as a thermal bath - can be used to correctly de- 
termine the "transition rates", i. e. that quantity that 
control the inter-basins dynamics. 

According to the second point of view, the expression 
of the free energy F{uj, ksT, ...) is found integrating over 
the degrees of freedom orthogonal to w, and the result 
is the usual free energy in the canonical ensemble. This 
free energy defines the pure states as local minima of the 
free energy (a pure state is now a not necessarily com- 
pact set of points in the configuration space); the free 
energy harriers among these pure states signal transitions 
among states. Yet, once a coarse-grained description of 
the system in terms of reaction coordinates uj is given, 
the possibility of discriminating among energy and en- 
tropy barriers is found. As a back-of-the-envelope calcu- 
lation to illustrate the concept, let's consider a collection 
of iV harmonic degrees of freedom plus a double well de- 
gree of freedom with two equal minima separated by an 
energy barrier of height A. The reaction coordinate cor- 
responds obviously to the anharmonic degree, and the 
transition rate among the minima can be calculated to 
be oc exp{—A/kBT). If, on the contrary, we work in the 
spirit of the second point of view, i. e. with trajectories 
at constant potential energy U, we can estimate the tran- 
sition rate as a ratio of volumes, and specifically as the 
ratio of the volume ~ (Vu)^ relative to the time that 
the system spend in the basin of a minimum and the 

volume ~ (^^/{U ~ A)j pertaining to the time spent 

in the order-one saddle basin. In the thermodynamic 
limit, where N ^ oo and U NkBT/2, we have that 
the transition probability (1 — iA/NksT)^^'^ converges 
to exp(—A/kBT). We can conclude that the two points 
of view give rise to the same expression for the transition 
probability. 

As a further example we present the numerical com- 
parison between the " true" dynamics of a finite size sys- 
tems and a fictions dynamics where the same system is 
forced to evolve at constant potential energy. In particu- 
lar, we studiy the self diffusion coefficient in the so called 
" Trigonometric Model" (TM) , introduced by Madan and 
Kcyes as a dynamical system with a very simple 
structure of the PES. In the TM each degree of freedom 
is independent from the other and experiences a periodic 
external force, in particular the TM potential energy has 
the form: 

N 

Vi{x}^=l..N) = ^ A[l - COs{TTX^/do)]. (5) 
i=l 

Despite its simplicity, the PES of this model contains 
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some features that have been retrieved in more reahstic 
glasses. The TM shares with the binary mixture - and 
with the modified - LJ system investigated in Ref. the 
existence of a regular organization of the saddles above 
a given minimum (the elevation in energy of the sad- 
dle is proportional to the saddle's order) and a regular 
distribution of the minima in the configurational space 
(nearest neighbor minima lie at a well defined Euclidean 
distance). What the TM misses is the existence of min- 
ima with different energies, a characteristic of realistic 
glasses [Q. This deficiency is not important here, as we 
want only to select a toy model to compare the "canon- 
ical" and " iso-potential" dynamics. In order to simplify 
the calculation, in this example we treat the dynamical 
evolution of the system following a Langevin dynamics 
in the high friction limit. 

In a first step we evaluate the time evolution of the N 
variables Xi{t) according to the true Langevin dynamics: 



(6) 



where r]i{t) is a random variable, 5-correlated in time 
and with variance proportional to the temperature, < 
Vi{t)'njit') >= '2S^jS{t - t')TkBT, and T fixes the time- 
scale of the evolution. In the case of the stochastic motion 
in periodic potential (as the TM), the self diffusion coef- 



ficient D = limt^oo{N t) 
be calculated analytically 



\xi{t) - Xi{Q)Y > can 



and turns out to be: 
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r 1 



/3 Z(/3) Z(-/3) 



(7) 



where Z{/3) is the partition function for a single degree of 
freedom. For potentials - like that of Eq. (||) - having a 
skew symmetry, V{x) = 2A — V{x — do/2), equation (|^) 
simplifies to: 



(8) 



Finally, by calculating the partition function for the TM, 
one gets 



(9) 



where Io{x) is the Bessel function of order 0. The dimen- 
sionless quantity D/Dfree, where Dfree — ^ksT is the 
diffusion coefficient for free brownian motion, is reported 
as a function oi 9 = kuT j A in Fig. ^ (full line). 

As second step, we calculate the diffusion coefficient 
for the same potential model and the same stochastic 
dynamics, but with a constrain that obliges the system 
to move at constant potential energy. The motion at 
constant potential energy implies that iiV^V" — 0, 
and this condition is automatically satisfied by modifying 
the equation of motion (0) as: 



Xi(t) 



(10) 



The evolution of iV=1000 variables obeying Eq. ( pX)| ) has 
been solved numerically at different temperatures. The 
diffusion coefficients calculated in these runs (for each T 
we made 500 replica each one 500' 000 time steps long) 
are reported in Fig. ^ (open dots). The comparison be- 
tween the values of D evaluated along canonical tra- 
jectories (full lines) and along iso-potential trajectories 
(open circles) clearly indicates that, already for not-too- 
large system sizes (A^=1000), the "true" dynamics take 
place along constant potential energy manifolds, and that 
the potential energy fluctuations, that are present in the 
canonical calculation, are not relevant for describing the 
diffusion process. 

As we have shown, and substantiated with the previ- 
ous examples, in the thermodynamic limit the dynam- 
ics of the system as a whole takes place along constant 
potential energy manifolds, the fluctuation becoming ir- 
relevant. In this framework, i. e. discussing the under- 
cooled dynamics as motion of the representative point 
in the configuration space, the notion of activated pro- 
cess becomes irrelevant as well. Either at high and low 
temperature, above or below Tmcti the representative 
point in the configuration space never crosses potential 
energy barriers. The slowing down of the dynamics with 
decreasing T is a consequence of a topological properties 
of the PES, i. c. the decreasing number of the accessible 
iso-potential paths between minima with decreasing po- 
tential energy. Some of these open paths still exist even 
on the typical constant potential energy manifolds visited 
at temperature well below Tmct, as shown, for example, 
in Ref. ||ll[] . One can conclude that the dynamical slowing 
down is fully entropy controlled and no qualitative differ- 
ences in the way the system changes basins are observed 
crossing Tmct- On the contrary, a qualitative change in 
the dynamic is expected to take place at even lower T (at 
the Kautzmann temperature Tk), where the number of 
accessible basins becomes non extensive. In this respect, 
whether or not Tk is different from zero is still matter of 
debate. 
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FIG. 1. a) Sketch of the potential energy surface (PES) 
along a direction with negative curvature that could lead to 
inter-basin diffusion. The marked point represents the posi- 
tion of the system; for T > Tmct the apparent energy bar- 
rier between this point and the top of the PES along the 
direction considered is less than ~ 10~^ ksT for the models 
studied. Although the section of the PES in this example has 
two minima along the direction of the eigenmode associated 
with this negatively curved direction, it can be that these ap- 
parent minima lead to the same basin b) The fraction 
< fdiff > of diffusive modes in the SPC/E model for wa- 
ter goes to zero at the mode-coupling transition temperature 
Tmct (from ref. M). 
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FIG. 2. Effect on the landscape from switching from the 
potential U to the pseudo-potential |VC/p: the saddle S be- 
comes a minimum with a finite basin of attraction, the basins 
of attraction of the minima A and B become smaller. Note 
that the picture is a two-dimensional sketch, so that the bor- 
der of the basins "acquires" a small volume: in higher dimen- 
sions the effect is enanced. 
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FIG. 3. a) Temperature dependence of the potential en- 
ergy of the system, of the saddles and of the minima for the 
Modified Lennard- Jones system of ref. The full squares 

indicate extrapolations of the energies at T = Tmct- Note 
that, in the range of the graph, U(Tmct) is still much above 
the potential energy range of the saddles, b) Schematic repre- 
sentation of the motion of the system at a temperature slightly 
below Tmct- For a big system, the motion happens essen- 
tially along constant potential energy. Although there are 
still saddles below the trajectory, the system is seldom near 
the border above the saddles; when the trajectory meets the 
border, a change of basin happens without energy activation. 
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— canonical ensemble (theoretical) 
O constant potential energy (numerical) 



FIG. 4. Diffusion in the trigonometric model The 
continuous hne is the analytical solution for the dynamics in 
the canonical ensemble ; circles are numerical calculations 
for the constant potential energy dynamics (Eq. 
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